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We introduce a major theoretical generalization of existing techniques for handling the three- 
body problem that accurately describes the interactions among four fermionic atoms. Application 
to a two-component Fermi gas accurately determines dimer-dimer scattering parameters at finite 
energies and can give deeper insight into the corresponding many-body phenomena. To account for 
finite temperature effects, we calculate the energy-dependent complex dimer-dimer scattering length, 
which includes contributions from elastic and inelastic collisions. Our results indicate that strong 
finite-energy effects and dimer dissociation are crucial for understanding the physics in the strongly 
interacting regime for typical experimental conditions. While our results for dimer-dimer relaxation 
are consistent with experiment, they confirm only partially a previously published theoretical result. 
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The physics of strongly interacting fermionic sys- 
tems is of fundamental importance in many areas of 
physics encompassing condensed matter physics, nuclear 
physics, particle physics and astrophysics. The last 
few years have seen extensive theoretical and experi- 
mental efforts devoted to the field of ultracold atomic 
Fermi gases. The ability to control interatomic in- 
teractions through magnetically tunable Feshbach res- 
onances has opened up broad vistas of experimentally 
accessible phenomena, providing a quantum playground 
for studying the strongly interacting regime. For in- 
stance, near a Feshbach resonance between two dissimilar 
fermions, the s-wave scattering length a can assume pos- 
itive and negative values, allowing for the systematic ex- 
ploration of Bose-Einstein condensation (BEC) and the 
Bardeen-Cooper-Schrieffer (BCS) crossover regime, in 
which bosonic (a > 0) and fermionic (a < 0) types of su- 
perfluidity connect smoothly [1, 2]. In this broad context, 
few-body correlations [3, 4] play an important role in de- 
scribing the dynamics of such systems. On the BEC side 
of the resonance (a > 0), dissimilar fermions pair-up into 
weakly-bound bosonic dimers, and the zero (collision) en- 
ergy dimer-dimer scattering length, a dd (0), determines 
various experimental observables such as the molecular 
gas collective modes, the internal energy, and even the 
macroscopic spatial extent of the confined cloud [1, 2]. 
Although a better description of the many-body behav- 
ior has emerged through the inclusion of few-body corre- 
lations, most of the current understanding of crossover 
physics relies on zero-energy theories, and very little 
is known about finite energy effects in this regime (see 
Ref. [5] and references therein). 

In this Letter we demonstrate important finite energy 
effects which can potentially impact the physics of a finite 
temperature ultracold Fermi gas in the crossover regime. 
Our results show deviations from zero-energy dimer- 
dimer collisions and indicate that, at experimentally rele- 
vant temperatures and scattering lengths, molecular dis- 
sociation might play an important role. The crossover 
regime can be viewed as a long-lived atom-molecule mix- 
ture, where dimers are dynamically converted to atoms 
and vice-versa. In order to account for finite tempera- 



ture effects, we calculate the energy dependent complex 
dimer-dimer scattering length, add(E co i), where E co \ is 
the collision energy. The real and imaginary parts of 
add correspond, respectively, to contributions from elas- 
tic and inelastic (dissociative) collisions [6], both of which 
should be considered to properly model the Fermi gas at 
realistic temperatures. In the zero-energy limit we repro- 
duce the well known prediction ow(0) ~ 0.6a [4]. How- 
ever, when the dimer binding energy, Ef, = ft 2 /(2/i2b& 2 ) 
(where /i 2 b is the two-body reduced mass) is compara- 
ble to the gas temperature T, finite energy effects and 
molecular dissociation become important, defining a crit- 
ical scattering length a c = ft/(2/Z2b&B7 1 ) 1//2 , where ks is 
Boltzmann's constant, beyond which an atom-molecule 
mixture should prevail. 

We also study dimer-dimer relaxation, in which two 
weakly-bound dimers collide and make an inelastic tran- 
sition to a lower energy state. In such a process, the 
kinetic energy released is enough for the collision part- 
ners to escape from typical traps. Ref. [4] predicted 
that, near a Feshbach resonance, dimer-dimer relaxation 
is suppressed as a -2 ' 55 , explaining the long lifetimes ob- 
served in several experiments [2, 7]. Here we also verify 
this suppression, although with an a dependence that is 
not described as a simple power-law scaling as originally 
predicted [4]. While the a -2 55 scaling law has already 
been tested (Regal et al. [7] found a - 2 - 3 ±°- 4 and Bourdel 
et al. [2] a~ 2 - 0±0 - 8 ), our calculations demonstrate that 
finite range corrections can explain the apparent experi- 
mental scaling law behavior, despite deviations from that 
power-law for larger a. 

We solve the four-body Schrodinger equation in the 
hyperspherical adiabatic representation, which offers a 
simple yet quantitative picture. A finite range model is 
assumed for the interatomic interaction, and a physically- 
motivated variational basis set is adopted to solve the hy- 
perangular equations [8] . While several hyperangular pa- 
rameterizations exist, we find that the best choice is the 
"democratic" hyperspherical coordinates [9] in which all 
possible fragmentation channels are treated on an equal 
footing, which describes elastic and inelastic processes in 
an unified picture. 
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In the adiabatic hyperspherical representation, the col- 
lective motion of the four fermions is described in terms 
of the hyperradius R, characterizing the overall size of the 
system. The interparticle relative motion is described by 
the hyperangles = {6*i, 9 2 , (j>i, (^3}, and the set of 
Eulcr angles {a, ft, 7} specifying the orientation of the 
body-fixed frame [9]. Q\ and 62 parameterize the mo- 
ments of inertia while <j)\ , 4>2 and parameterize inter- 
nal configurations [9]. Integrating out the hyperangular 
degrees of freedom, the Schrodinger equation reduces to 
a system of coupled ordinary differential equations, given 
in atomic units (used throughout this Letter) by: 
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F u (R)+Y,W^(R)F u ,(R) = 0, (1) 



identical fermions and assume the one between dissimilar 
fermions to be v(r) — Dsech 2 (r /ro) , where r is the inter- 
atomic distance and D is tuned to produce the desired 
a. We choose the atomic mass m and effective range 
r = 181 a.u. [11] to be those of 40 K. The eigenval- 
ues and eigenfunctions of -ff a d, namely the hyperspher- 
ical potentials U„(R) and channel functions <I>„(i?;f2), 
determine the effective potentials and nonadiabatic cou- 
plings in Eq. (1): W vv = U v - Q vl/ j2[i and W vv > = 
-[P vu >d/dR+Q^]/2fi, where P vv , = ($ u \d/dR\$„>) 
and Q vv > = {<S> v \d 2 /dR 2 \<S> v .). We find <S> v {R]ti) varia- 
tionally by expanding in exact eigenfunctions of Eq. (2) 
at large and small R [8]. At ultracold energies the con- 
vergence of the scattering observables with respect to the 
number of basis functions is surprisingly fast [8] . 



where = to/43 is the four-body reduced mass (to being 
the atomic mass) , E is the total energy, F u is the hyperra- 
dial wavefunction, and v represents all quantum numbers 
needed to label each channel. Scattering observables can 
then be extracted by solving Eq. (1), where the nonadia- 
batic couplings W u:j L u i drive inelastic transitions between 
channels described by the effective potentials W vv . 




FIG. 1: (color online). Analysis of the probability density 
integrated over the hyperangles d\ and 62 at R = 0.41 a 
is shown, (a) An isosurface of the probability density at 
j <E> (7?; O) 1 2 = 0.1|$(.R; f2)|max f° r the dimer-dimer channel is 
shown. The darker (lighter) colors correspond to a more (less) 
linear configuration for the four-particle system. Q(x) is the 
unit-step function, (b)-(e) show density plots for fixed values 
of 03. The darker regions represent higher probabilities for 
which planar configurations are shown to illustrate the most 
probable four-body geometry at selected points. 

In the hyperspherical representation, the major reduc- 
tion to Eq. (1) is accomplished by finding eigenfunctions 
of the (fixed R) adiabatic Hamiltonian, 



H ad (R,n) 



k 2 (n) + i2 

2[iR 2 



v(R,n). 



(2) 



In the above equation, A is the grand angular momen- 
tum operator [9] and V includes all two-body interactions 
[10]. For simplicity, we neglect the interaction between 
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FIG. 2: (color online), (a) Several four-fermion hyperspheri- 
cal potentials attached to all relevant breakup thresholds are 
shown, (b) The energy-dependent elastic (red) and inelastic 
(green) parts of add [Eq. (3)] are shown. For energies £<£i 
we find a r dd = 0.605(5)o [4], while for E co i m Et finite energy 
corrections strongly affect add- Solid black line: add obtained 
from the effective range expansion [12]. 

However, including higher order correlations that de- 
scribe dimer-atom-atom and four free atom configura- 
tions is crucial for accurately describing scattering pro- 
cesses at any collision energy. We find that for R< a the 
strongest contribution to the probability density of the 
dimer-dimer channel function comes from dimer-atom- 
atom like configurations. Figure 1 shows a graphical rep- 
resentation of this channel function in terms of the inter- 
nal configuration angles <f>\, 4>2 and 4>3- The four "cobra" - 
like surfaces explicitly illustrate the four-fold symmetry 
(52 <8> 52) of the fermionic problem. The "spines" of 
the cobras correspond to the interaction valleys where 
two dissimilar fermions are in close proximity while the 
"hoods" loosely represent the larger phase-space explored 
by dimer-atom-atom like configurations. 

Figure 2(a) shows the hyperspherical potentials for a = 
125ro showing the full energy landscape with the four- 
body thresholds for dimer-dimer (FF' + FF'), dimer- 
atom-atom (FF'+F+F'), and four atom (F+F'+F+F') 
collisions. Notice that the four-body potential associated 
with dimer-dimer collisions is repulsive for R < a, in- 
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FIG. 3: (color online), (a) The scattering length dependence of a r dd and a dd (inset) at fixed collision energies is shown. At any 
non-zero collision energy a dd deviates from the zero-energy prediction (black dashed line) as a — > oo and strong contributions 
from molecular dissociation, a dd , occur, (b) The dimer-dimer relaxation rate, V^f , is shown versus a. The solid line is the total 
V^f and the dashed lines are contributions from different decay pathways (see text). For intermediate a, we reproduce both 
experimental data [7] (filled circles) and the a -2 ' 55 scaling law [4] while deviating from that for larger a. Inset: we show T p as 
a function of R for a = lOOro, 80ro, 65ro, and 50ro (red, green, blue, and magenta respectively). 



dicating that zero-energy dimer-dimer elastic scattering 
must be qualitatively similar to scattering by a hard- 
sphere of radius a, i.e., a dd (0) oc a. Although a clear and 
qualitative picture emerges from the four-body potentials 
alone, we in practice extract scattering observables from 
coupled-channel solutions to Eq. (1). We define the en- 
ergy dependent dimer-dimer scattering length, a^(E c0 {), 
in terms of the complex phase-shift obtained from the 
corresponding S- matrix element [Sdd,dd = exp(2i^d)] 5 

add(E C oi) = : — — = a r dd (E co i) + i a l dd {E col ). (3) 

Here, k dd = 2mE co i, E co i = E + 2Eb is the collision en- 
ergy and a dd and a dd < are the real and imaginary parts 
of a dd , representing clastic and inelastic contributions [6] . 

Figure 2(b) shows a r dd and a dd for a = 125r . For 
energies E co i -C £), we find that add(0) = 0.605(5)a in 
agreement with Refs. [4, 12] while for E co i < al- 
though molecular dissociation is still not allowed, i.e. 
a dd ~ 0, we obtain strong corrections to the zero-energy 
result. At these energies, an effective range expansion, 
a dd( E coi) = Odd(0) - \r d dk 2 dd where r dd = 0.13a [12], 
is accurate over a small range, but quickly fails to re- 
produce our results [see black solid line in Fig. 2(b)]. 
For E co i > Eb, the channels for molecular dissociation 
become open leading to strong inelastic contributions 
to add(Ecoi), as parametrized by a dd . Our results in- 
dicate that both a dd and a dd are universal functions of 
energy and scattering length, i.e., insensitive to the de- 
tails of the short-range physics, which should extend up 
to E co i <C 1/mr-g in the absence of deeply bound states. 
Due to the small number of basis functions used in these 
calculations, our results for Eb -C E co i <C l/(mrg) are 
not fully converged, but we expect their qualitative be- 
havior, i.e., the sharp decease in add{E co i), to persist. 



Figure 3(a) demonstrates that when approaching the 
Feshbach resonance (a— >oo) at any finite collision energy, 
molecular dissociation becomes increasingly more impor- 
tant and a r dd substantially deviates from the zero-energy 
predictions [black dashed line and inset in Fig. 3(a)]. As 
a — > oo, Eb cx 1/a 2 becomes extremely small and such 
finite energy effects [see Fig. 2(b)] are relevant even at ul- 
tracold energies. Therefore, the molecular binding energy 
Eb, or equivalently a c = 1/ ^/Tfl^T, defines the range be- 
yond which (i.e., a > ac ) deviations from the zero-energy 
predictions can be observed. Perhaps more importantly, 
it specifies a regime beyond which molecular dissociation 
can lead to a long-lived atom-molecule mixture [13, 14], 
where dimers are continuously converted to atoms and 
vice- versa, i.e. FF' + FF' <h> FF' + F + F' . Further, 
this indicates that the underlying physics of the strongly 
interacting regime may fundamentally depend on tem- 
perature. Values for a c at 100 nK are 7000 a.u. for 40 K 
and 17000 a.u. for 6 Li, and therefore the finite energy 
effects above can become experimentally relevant [2]. 

We also study vibrational relaxation due to dimer- 
dimer collisions. We verify the suppression of the re- 
laxation rate as a — > oo, however, with a different a de- 
pendence than a~ 2,55 predicted in Ref. [4]. In Ref. [4] 
it was assumed that the main decay pathway for relax- 
ations is a purely three-body process and requires only 
three atoms to be enclosed at short distances. Therefore, 
it neglects the effects of the interaction with the fourth 
atom. Here, however, we analyze such effects and find 
that it strongly influences the suppression of relaxation. 
In our calculations we express the inelastic transitions 
probability T p (a, R) in terms of the probability of having 
three atoms at short distances as a function of the dis- 
tance m R of the fourth atom from the collision center 
[15]. We calculate T p from our fully coupled-channel so- 
lutions and effective potentials [see Figs. 1 and 2(a)] and 
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our results are shown in the inset of Fig. 3(b). 

In our model the relaxation rate is simply proportional 
to the transition probability T p (a,R). It is interesting 
note that our formulation allows for the analysis of dif- 
ferent decay pathways. For instance, at short distances, 
R w ro, T p describes inelastic transitions in which all 
four atoms are involved in the collision process. At large 
distances, R/a » 1, T p describes the decay pathway 
where only three atoms participate in the collision, akin 
to the process studied in Ref. [4]. We note, however, 
that for values of R up to R m 5 the scaling law for 
relaxation depends strongly on R/a and greatly deviate 
from the a~ 2 ' 55 scaling. In order to take into account 
inelastic processes for all values of R we define an effec- 
tive transition probability by integrating T p (a, R) over R 
[15]. Our results for the relaxation rate, V^f, are shown 
in Fig. 3(b) where the red solid line is obtained by in- 
tegrating T p from R = 2r up to 10a [16] giving an ap- 
parent scaling law of a - 3 - 20 ± 05 . The dashed lines are 
obtained from integrating T p from R — 2ro to 5ro and 
from R = 5r to 10a, which yields scaling laws of a -4 02 
and a -3- 20 ±0- 05 ; respectively, "separating" the contribu- 
tions from the decay pathways in which four and three 
atoms participate in the collision process. The ampli- 
tudes for each of these contributions, however, are dis- 
connected as they depend on the details of the four- and 
three-body short-range physics. In contrast, the ampli- 
tudes for the a~ 3,20 and a -2 55 processes are governed by 
the same three-body physics. As a result, the fact that 
we don't observe the a -2 55 scaling implies that it is not 
important for the range of a used here. The amplitude 
for the process which leads to the a -2 55 scaling is expo- 
nentially suppressed owing to the unfavorable overlap of 



the dimers' wavefunction [see inset of Fig. 3(b)]. In fact, 
for our largest values of a, it is already apparent that 
in the very large a limit the rate deviates from a -3 ' 20 , 
however, to a behavior different than a~ 2,55 [15]. 

Figure 3(b) rescales our results for V^f by an overall 
constant chosen to fit the experimental data for 40 K at 
a temperature of 70 nK (Regal et al. [7]). We note, 
however, that between a = 1000 and 3000 a.u. [17], our 
results agree with both the experimental data and the 
fl -2.55 sca ii n g i aW; approaching our predicted scaling law 
fl -3.20 on jy f or i ar g er values of a. This change in behavior 
of V^f originates in the finite range of our model, which 
represents physics beyond the zero-range model of Ref. [4] 
where the a~ 2 ' 55 scaling applies for all a. 

In summary, we have calculated the energy depen- 
dent complex dimer-dimer scattering length, add(E co i), 
by solving the four-body Schrodingcr equation in the adi- 
abatic hyperspherical representation. Our results demon- 
strate that for experimentally relevant temperatures and 
scattering lengths the clastic and inelastic contributions 
of add are equally important. We show that molecular 
dissociation plays an important role and suggest that the 
many-body behavior in the strongly interacting regime 
might be significantly altered at finite temperature. Our 
results also demonstrate a stronger suppression for dimer- 
dimer relaxation, compared to that obtained in Ref. [4], 
while remaining consistent with experimental data. 
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ence Foundation. 



[1] D. M. Eagles, Phys. Rev. 186, 456 (1969); A. J. Leggett, 
J. Phys. C 41, 7 (1980); M. Holland et al, Phys. Rev. 
Lett. 87, 120406 (2001); E. Timmermans et al., Phys. 
Lett. A 285, 228 (2001); Y. Ohashi and A. Griffin, Phys. 
Rev. Lett. 89, 130402 (2002). 

[2] C. Chin et al, Science 305, 1128 (2004); T. Bourdel et 
al, Phys. Rev. Lett. 93, 050401 (2004); M. Bartenstein 
et al, ibid. 92, 120401 (2004); C. A. Regal, M. Greiner, 
and D. S. Jin, ibid. 92, 040403 (2004); M. W. Zwierlein 
et al, ibid. 92, 120403 (2004); G. B. Partridge et al, 
ibid. 95 020404 (2005); M. Bartenstein et al., ibid. 92, 
203201 (2004); J. Kinast et al, ibid. 92, 150402 (2004). 
M. Greiner, C. A. Regal, and D. S. Jin, ibid. 94, 070403 
(2005); M. W. Zwierlein, et al, Nature (London) 435, 
1047 (2005). 

[3] P. Pieri and G. C. Strinati, Phys. Rev. B 61, 15370 
(2000). 

[4] D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. 
Rev. Lett. 93, 090404 (2004); Phys. Rev. A 71, 012708 
(2005). 

[5] M. J. Wright et al, Phys. Rev. Lett. 99, 150403 (2007). 
[6] J. L. Bohn and P. S. Julienne, Phys. Rev. A 56, 1486 
(1997); R. C. Forrcy et al, ibid. 58, R2645 (1998). 



[7] J. Cubizolles et al, Phys. Rev. Lett. 91, 240401 (2003); 

K. E. Strecker, G. B. Partridge, and R. G. Hulet, ibid. 

91 080406 (2003); C. A. Regal, M. Greiner, and D. S. 

Jin, ibid. 92, 083201 (2004). 
[8] Nirav P. Mehta, Seth T. Rittenhouse, Jose P. D'Incao, 

Chris H. Greene, arXiv:0706.1296. 
[9] V. Aquilanti and S. Cavalli. J. Chem. Soc. : Faraday 

Transactions 95, 801 (1997); A. Kuppcrmann. J. Phys. 

Chem. A 101, 6368 (1997). 
[10] We omit all reference to the Euler angles because this 

study treats only the dominant partial wave J 7 * — + , 

where J is the total orbital angular momentum and tt is 

the parity quantum number. 
[11] V. V. Flambaum, G. F. Gribakin, and C. Harabati, Phys. 

Rev. A 59, 1998 (1999). The value for r is obtained 

from the expression r = 2 1/2 3" 1 r(i)r(|)" 1 (mC 6 ) 1/4 , 

Eq, (24), at resonance. 
[12] J. von Stecher, Chris H. Greene, and D. Blume, Phys. 

Rev. A 76, 053613 (2007); ibid. 77, 043619 (2008). a dd 

in these works is calculated to an accuracy similar to our 

present value. 

[13] S. Jochim et al, Phys. Rev. Lett. 91, 240402 (2003). 
[14] Y. Shin, A. Schirotzek, C. H. Schunck, and W. Ketterle, 



Phys. Rev. Lett. 101, 070404 (2008). 
[15] See Supplemental Material attached to this file for a more 
detailed description of our model for dimer-dimer relax- 
ation. 

[16] We have found that integrating T p from R = 2ro up to 
10a is enough to ensure that contributions for R < 2r 
and R > 10a are neglegible. 

[17] The experimental data for a > 3000 a.u. were taken in 
the regime where the molecular size is expected to be 
larger than the interparticle spacing, which prohibits a 
proper comparisson to our results. 



Supplemental Material: "Dimer-dimer collisions at finite energies in two-component 

Fermi gases" 

J. P. D'Incao, Seth T. Rittenhouse, N. P. Mehta, and Chris H. Greene 
Department of Physics and JILA, University of Colorado, Boulder, Colorado 80309-0440, USA 



In this supplemental material we provide additional 
details of our model for dimer-dimer relaxation in terms 
of the transition probability obtained from our numerical 
calculations [1]. 

The key observation in our model is that the inelastic 
transitions leading to deeply bound molecular final states 
can only occur when at least three atoms are enclosed at 
distances comparable to tq. In the hyperspherical repre- 
sentation, the decay pathway in which only three atoms 
participate is viewed as an infinite series of avoid crossing 
between the initial dimer-dimer channel and all possible 
final states, as illustrated in Fig. 1(a), where the red solid 
curve describe the initial collision channel and the green 
dashed curves some of the possible final states. 

Therefore, for a given a, the inelastic transition to a 
particular final state can be described in terms of the 
Fermi Golden rule 

rW(o) ex \{* M (r, si) \v (R, n)\ * A (r, n))\ 2 (i) 

where *f?\ is the final state wave function, labeled by the 
quantum number A, ^ dd is our fully coupled dimer-dimer 
wave function, and V is the sum of the interatomic inter- 
actions. The hyperangular behavior of the integrand is 
assumed to be proportional to the probability amplitude 
of three particles being in close proximity. On the other 
hand, the hyperradial behavior of the outgoing channel 
will oscillate very quickly away from the classical turning 
point, as shown schematically in Fig. 1(b). The fast os- 
cillation will, in general, cancel out in regions away from 
the classical turning point, R\. This indicates that the 
integral over the hyperradius will be proportional to the 
area within the first oscillation of ^\ (R) times the re- 
maining hyperradial behavior evaluated at R\ . Based on 
these considerations, the squared result yields 
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T^\a,R x ) <x -^\F dd (Rx)\ z F(Ri), 



(2) 



where Fdd is the dimer-dimer hyperradial wave function 
and T (R\) is the probability of having three out of four 
atoms at distances comparable to ro at hyperradius R\. 
Therefore, the inelastic transitions to a particular final 
state A occurs in the vicinity of R\ and it is governed 
by the probability of having the three atoms within dis- 
tances comparable to ro. 

In practice, we have calculate J 7 (R\) by defining the 
proximity operator 



f(R,Q) = e 



- p -(rt 2 +rL)/2r 2 



cyclic permutaions, (3) 
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FIG. 1: (color online), (a) A schematic representation of 
the dimer-dimer relaxation process is shown. The red-solid 
solid curve represents the effective potential for the initial 
dimer-dimer channel, FF' + FF', and green-dashed curves 
represent the possible final decay channels involving at least 
one deeply bound molecular state FF . Inelastic transitions 
to this almost continuum of final states are allowed for all 
values of R (see main text), (b) Qualitative representations of 
the behavior of the hyperradial solutions in the dimer-dimer 
channel (black dot-dashed curve) and in a particular final 
channel (solid black curve) are shown to illustrate that the 
inelastic transition is more likely near the classical turning 
point R\. 



close to each other. V is then simply defined as 



T{R) a (<l>dd(R;n)\f(R,n)\<S>dd(R;n)), 



(4) 



which is non-zero only when three atoms are sufficiently 



where $dd is our fully coupled dimer-dimer channel func- 
tion and the integration is taken over all the hyperangles. 

Our model for dimer-dimer relaxation, therefore, is 
simply obtained by summing Eq. (2) over the near con- 
tinuum of A states, approximated by an integral over the 
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classical turning points, which yields 
Kef « ^- / Tf>{a,R x )p(R)dR 

= ^ [ ^\F dd (Rx)\ 2 T(Rx) P (Rx)dR x 
Kdd J H\ 

_ 1 / PwKB(Rx)f(Rx) fry \jt> CO 

~ i— / ^ — 7^ piRxjdRx, (5) 

fcdd J RxK(Rx) 

where fc 2 d = 2m(E + 2£'f ) ), p (-Ra) is the nearly constant 
density of states and \Fdd (Rx) | 2 was approximated by 
the WKB wavefunction in the classically forbidden re- 
gion, \F dd (Rx)\ 2 ~ Pwkb(Rx) /k(R x ). Here, P W kb 
is the WKB tunneling probability for the dimer-dimer 
hyperradial wavefunction and k wkb (Rx) is the WKB 
wave number. 




a (a.u.) 



In our model for relaxation, breaking up the integra- 
tion over T p (a) [Eq. (5)] allows us to analyze the contri- 
butions from different pathways. In Fig. 2 we show the 
total rate as a solid-red curve, obtained by integrating 
Eq. (5) from 2tq to 10a [2]. In addition to that, we also 
plot the results for V^f obtained by integrating Eq. (5) 
from 2ro to 5ro, see green-dashed curve in Fig. 2. This 
result determines the contribution from inelastic tran- 
sitions which occur predominantly when all four atoms 
are within distances comparable to r . The blue-dashed 
curve in Fig. 2, however, shows our results obtained by 
integrating from a to 4a, determining the behavior of 
the contributions from inelastic transitions that occurs 
near R — a. The main difference between this contri- 
bution and the total rate comes from the inelastic tran- 
sitions for R < a. Although numerically we are unable 
to go to larger values of a, it is clear that the contri- 
butions for transitions near R — a becomes increasingly 
more important and in the very large a limit we expect 
these contributions to dominate the total rate. Interest- 
ingly, for the values of a we studied, the contribution 
for transitions near R — a already falls off slower than 
the a -2 ' 55 prediction of Ref. [3]. Therefore, we conclude 
that the mechanism that leads to the a -2 55 suppression, 
although, qualitatively correct, doesn't quantitatively de- 
scribe the transitions near R = a due to the presence of 
the fourth atom. 



FIG. 2: (color online). The vibrational relaxation rate for 
different decay pathways is shown. The red-solid curve is 
the total rate given by Eq. (5), the green-dashed curve is 
the contribution from short-range inelastic transitions where 
all four-atoms participate in the collision. The blue-dashed 
curve is the contribution from inelastic transitions near R = a 
involving only three-atoms illustrating the effects due to the 
presence of the forth atom (see main text). 
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